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We investigate the constraints on the scalar, vector and spin-3/2 dark matter interaction with the 
standard model particles, from the observations of dark matter relic density, the direct detection 
experiments of CDMS and XENON, and the indirect detection of the p/p ratio by PAMELA. A 
model independent way is adopted by constructing general 4-particle operators up to dimension 
... 6 for the effective interaction between dark matter and standard model particles. We find that 

the constraints from different experiments are complementary with each other. Comparison among 
, these constraints may exclude some effective models of dark matter and limit some parameters of 

. others. The spin-independent direct detection gives strong constraints for some operators, while the 

indirect detection of p/p data can be more sensitive than direct detection or relic density for light 
I dark matter (whose mass < 70 GeV) in some cases. The constraints on some operators for spin-3/2 

■ dark matter are shown to be similar to those on their analogous operators for Dirac fermionic dark 

matter. There are still some operators not sensitive to the current dark matter direct and indirect 
I search experiments. 
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I. INTRODUCTION 



' It is by now established [ll-Q that about 23% of the constituents of the Universe are composed of dark matter 
• (DM) [Jill- However, its nature remains unclear. A well- motivated candidate for DM is the weakly interacting 
(N| ' massive particle (WIMP), which must be stable, nonrelativistic, electrically neutral and colorless. If the WIMP mass 
lO . is from a few GeV to TcV while their interaction strength is of the weak scale, they can naturally yield the observed 
O ■ relic density of DM [1]. Ahhough there is no stable WIMP in the Standard Model (SM), WIMP candidates exist 
, in various theoretical models tryin g to solve the SM problems at the weak scale, such as supcrsymmetric models 
CN ' Isl-fisj , extra dimensional models |1^ - [2H , little Higgs models [22| - [26j , left-right symmetric models pvl - lsoj , and some 
other models (e.g. [Hl-llll). 

The specific models mentioned above are very attractive, but still lack experimental support. The well-running 
LHC experiment may find some important signals of these models in the near future. However, if other new particle 
species are all so heavy that the DM particle is the only new particle within the reach of LHC, it will be very difficult 
to know which model the DM particle belongs to. In addition, it is possible that the DM particle may be first 
observed by direct or indirect detection experiments. These early observations may only provide information about 



' some general properties of the DM particle, and may not be able to distinguish the underlying theories. Therefore, the 
model-independent studies of the DM phenomenology can play an important role as they may avoid theoretical bias 
[ssl^ioj . Recently there have been quite a few papers to study various phenomenologies related with DM in the model- 
independent way l4l|-{58j . Especially the relic density measured by WMAP direct detection from CDMS Isol l60j 
and XENON (6ll. |62|. and possible collider signals from Tevatron (e.g. (63l. |64| ) and from LHC (e.g. [6^ |66| ) are 
considered in these studies. 

In our previous work (50| , we investigated a general set of 4-fermion operators for the effective interaction between 
the spin- 1/2 fermionic DM and the SM particles, and gave the phenomcnological constraints from the observed DM 
relic density, the direct detection experiments by CDMS and XENON, and the indirect detection of the p/p ratio 
by PAMELA [gJ. It was found that the constraints from different observations are quite complementary. Besides 
the possibility of spin-1/2 WIMPs, it is also possible that DM is composed of scalar, vector, or spin-3/2 WIMPs, 
which belong to different representations of Lorentz group. In these cases, the forms of possible effective operators are 
different. These differences may lead to distinguishable phenomcnological results. In this work we extend our previous 
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analysis to the cases of scalar, vector and spin-3/2 DM. Wc will consider a general set of 4-particle operators up to 
dimension 6 for the effective interaction between the WIMPs and the SM fermions and compute their phenomenological 
constraints. Wc will use the updated limit of XENONIOO SI direct detection [6l|, which is stronger than that adopted 
in our previous work jsof . 

This paper is organized as follows. In Sections HIl and Hill the effective models of scalar and vector DM are discussed, 
respectively. In the subsections of these two sections, we explore the constraints on these models from the DM relic 
density, direct and indirect detection searches and the validity of effective theory, and then present the combined 
constraints on the effective coupling constants of these models. In Section IIVI the study on the effective models of 
spin-3/2 DM are carried out briefly. The conclusions arc given in Section IVl 



II. SCALAR DARK MATTER 



Let us begin with the case that DM consists of complex scalar WIMPs (0 and its antiparticlc (f)^). We will add 
remarks when there is a notable difference between this case and the case of real scalar WIMPs. In order to study 
phenomenologies in a model-independent way, we construct effective interaction operators between the WIMPs and 
the SM particles. These interaction operators are limited only by the requirements of Hermiticity, Lorentz invariance 
and CPT invariance. 

We make the following assumptions similar to those in Ref. [H, [s^: (1) The WIMP is the only new particle species 
at the electroweak scale, and any other new particle species is much heavier than the WIMP. This implies that the 
thermal relic density of the WIMP is not affected by resonances or coannihilations. Thus it is possible to describe 
the interaction between the WIMPs and the SM particles in terms of an effective field theory. (2) The WIMP only 
interacts with the SM fermions through a 4-particle effective interaction of 4>''-4>-f-f type, but not with gauge or Higgs 
bosons. For simplicity, this interaction is assumed to be dominated by only one form in the set of 4-particle operators. 
(3) The WIMP-antiWIMP annihilation channels to the SM fermion-antifermion pairs dominate over other possible 
channels. In other words, the possible channels to final states that include gauge or Higgs bosons are assumed to be 
negligible. 

The effective interaction term of Lagrangian between two complex scalar WIMPs {(f) and (f)'') and two SM fermions 
(/ and /) is given by only one of the following expressions: 

Scalar int. (S) : = ^ 

5t0)/7V, (2) 



Vector int 


• (V): 




= E- 


Scalar-pseudoscalar int. 


(SP): 


^SP 


= E' 


Vector- axialvector int. 


(VA): 















^/2 



^4>^4>ni^L (3) 

^{ct^hX<^)h^^,f, (4) 



where the sum of / is over all the SM fermions, the symbol 9^ is short for 9^ — 9^, that is, i/)^ 9^ ^ = (9p(^) — {dfj,<f)^)(j), 
and the effective coupling constants Fj are real- valued numbers. Fj have mass dimension of —1 for the dimension-5 
operators (S and SP), and that of —2 for the dimension-6 operators (V and VA). We do not consider operators with 
more derivatives on the fields since these operators have higher mass dimensions such that they are more suppressed 
at low energy scale. Note that for the case of a real scalar WIMP, the WIMP vector current 0^c/) vanishes so 
that the V and VA effective operators vanish, leaving only the S and SP interactions. Since Ff are real-valued, the 
transformation properties of the 4-particle operators, Cs, Cy, Lsv and >CvA under C, P and T are the same as the 
corresponding 4-fermion operators of Ref. [501 ■ ^he operators are CPT invariant. 

Each form of the interaction operators listed above describes an effective model in which the WIMPs couple to 
the SM fermions. For each case, we will calculate the corresponding annihilation and scattering cross sections, which 
depend on the WIMP mass and the coupling constants Ff . Associated with the recent experimental results of the 
DM relic density, direct and indirect detection experiments, the phenomenological constraints on Ff can be derived, 
respectively. Since the sensitive regions of different kinds of experiments vary in different effective models, it would 
be interesting and meaningful to compare these constraints. 
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FIG. 1. The predicted coupling constants F[ as functions of the scalar WIMP mass M^, fixed by the observed relic density, 
Odm/i^ = 0.1109±0.0056 0], in each effective model. Frames in the left (right) column are results given for the case of universal 
couplings [Ff (X mj). 



A. Scalar WIMP annihilation and relic density 



In order to determine the relic density of WIMPs and the source function of cosmic-ray particles derived from 
WIMP annihilation in the Galactic halo, which is relevant to the DM indirect detection, we need to calculate the 
cross sections of WIMP-antiWIMP annihilation to fermion-antifermion pairs. For each case listed above, the result is 
given by 



CS.ann 



fV, ami 



fSP.ann 



fVA.ann 

where s is the square of the center-of-mass energy, is the WIMP mass, and c/ are the color factors, equal to 3 for 
quarks and 1 for leptons. For the case of real scalar WIMPs, due to the identity of WIMP and antiWIMP, the S and 
SP annihilation cross sections should include an additional factor of 4 if the Lagrangians are taken to be Eqs. ([T|) and 
([3]) with (j)'^ = (p. Eqs. ([S]) - ([5]) are also obtained in Ref. [s^. Taking the energy distributions of the initial particles 
into account, we need to calculate {(Jannv), the thermal average of cross section multiplied by relative velocity (or 
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FIG. 2. The predicted thermal reUc density (dashed Unes) of complex scalar WIMPs with S, V, SP and VA interactions 
respectively. In the left (right) column, results are given for the case of universal couplings (F/ cx m/). The narrow horizontal 
solid band shows the range of the observed DM relic density, Qdm/i^ = 0.1109 ± 0.0056 
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M0ller velocity, more exactly |68l|). 

In our previous work [sO], following Ref. we calculated (cTannw) in the "lab" frame (in which one of the two 
initial particles is at rest), because it can be proved that the calculation in this frame is equivalent to that obtained in 
the practical reference frame and more convenient. Thus we can expand s in the lab frame to be s ~ 4M| + M|w^ + 
jM^v"^ + 0{v^), and substitute it into Eqs. ([5]) - dS]). Then we have the form CTann^^ — Qo + o.i'v^ + 0{v'^), which leads 
to (aannv) — clq + 6aix^^ + 0{x^^) where x = M^/T and T is the temperature of DM. On the other hand, we can 
also directly calculate (aannv) in the practical reference frame following the method described in Ref. [6^ . These two 
methods lead to the same results, i.e.. 




(9) 
(10) 

(11) 
(12) 



Note that the above results disagree with those in Ref. [38 

The evolution of the DM abundance is described by the Boltzmann equation 



iHriA, = -((Ta 



cq CO 



(13) 



where H = a/a ^ ^J^irp/ (3Mpj) is the Hubble rate with Mpi denoting the Planck mass, {n^\) is the number 
density of WIMPs (antiWIMPs), and ('^^'t) ^h'^ corresponding equilibrium number density. For complex scalar 
WIMPs without particle-antiparticle asymmetry, we have = n^t . Thus the total DM particle number density is 
JT-DM = 2n^ [1, HI] . Using the standard procedure [1, Q to approximately solve the Boltzmann equation , we 
obtain a relic density of DM particles as 



DM' 



X 10-* GcV^ 



2.725 K 



Xf 



M.p\^/gJJf){aQ + iaix f ) 



(14) 



where Xf = M^/Tf with Tf being the freeze-out temperature, g*(Ty) is effectively relativistic degrees of freedom at 
freeze-out, Tq = 2.725 ±0.002 K [TO] is the present CMB temperature, and ao and ai are coefficients in the expansion 
CTann'i' — + fli + O (w** ) . The freczc-out temperature parameter Xf is evaluated numerically by solving the equation 



In 




45 gM^Mp\{ao + &aix 

Tj2 



2T:^^g^{xf)x^ 



(15) 



where the order one parameter c is taken to be 1/2 as usual. A notable difference between this calculation and our 
previous work [HO] is that the degree of freedom of a scalar WIMP or antiWIMP \s g = 1 , while that of a fermionic 
WIMP is g = 2. The numerical result of g*{T) is taken from Ref. [7l|. For the case of real scalar WIMPs, tidm = n^, 
and n^t in Eq. (jl3p should be replaced by so that the relic density will be one half of that for the complex scalar 
case given the same annihilation cross section. 

The 7-year observation of WMAP 0] gives the DM relic density nuufi^ = 0.1109 ± 0.0056. With this result we 
estimate the relation between the effective coupling constants Ff and the WIMP mass in each effective model, 
as shown in Fig. [TJ Two kinds of Ff are considered here. In the left column of Fig. [1] we show the results for the 
case when the effective couplings Ff to all the SM fermions are equal (universal couplings). In the right column of 
Fig. [TJ we show the results for the case of Ff ex m/ . This kind of proportionality may come from Yukawa couplings 
of a Higgs mediated interaction or some other unknown underlying mechanisms. In both cases, Ff decreases as 
increases for fixed JIdm^^ in each effective model. Besides, there are several interesting features in Fig. [TJ 

• In the case oi Ff (x m f , the curves of Ff vs. have obvious downward bends at about A/^ ^ rrit = 171 GeV. 



The reason is that when > nit , the annihilation channel 



tt is opened. Therefore because the couplings 
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are proportional to mj, this channel gives a tremendous contribution to the total (o'ann'^') 
gaps between t quark and other SM fermions. 



due to the huge mass 



• In the case of universal couplings, there are 2 pairs of nearly identical curves, that is, S ~ SP and V ~ VA if 
we denote this approximate identity of the two curves by a notation "~" for short here and henceforth. From 
Eq. ([9]) - ([T2|) we can see that in each pair their corresponding ((Tannf) differ only by terms of 0{T/M^) and/or 



terms of m'i/M^. These differences are not important when we calculate the cold relic density. 



• In the case of i^/ cx to/, S ~ SP and V ~ VA are still almost true except that in some small regions the deviations 
of the two nearly identical curves in each pair become larger. 

In Fig. [2] the curves of rt^Mh^ vs. for fixed coupling constants in the complex scalar WIMP models of S, V, 
SP and VA interaction operators are shown. In the left (right) column, results are given for the case of universal 
couplings {Ff (X TTif). The values of the couplings corresponding to the curves are denoted in each frames. All the 
curves in Fig.[5]bend more or less at about ~ 1.3, 4.2 and 171 GeV corresponding to the masses of c, b and t quarks 
respectively, due to the opens of these heavy quark annihilation channels. An important qualitative difference between 
the curves of effective operators with different mass dimension in both Figs. [T] and [2] is that the curves for dimension-6 
operators (V and VA) slope downward, while those for dimcnsion-5 operators (S and SP) remain horizontal in most 
mass range. The reason has been showed explicitly in Eqs. ([9]) - (fT2l). The annihilation cross section for dimension-6 
operators have additional factors of which suppress the relic density of large WIMP mass. 

It is worth noting that the results in Figs. [1] and [5] are based on the assumptions presented at the beginning of 
this section. Otherwise if resonances, coannihilations or annihilations to final states other than fermion-antifermion 
pairs are significant, the actual curves in Figs. [T] and [2] will be much lower than those shown there, as pointed out in 
Ref. M. 
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FIG. 3. The spin-independent (SI) WIMP-proton cross sections (daslied lines) for complex scalar WIMPs with S and V 
interactions. In the left (right) column, results are given for the case of universal couplings {Ff (x mj). The upper bounds set 
by CDMS II [5^ and XENONIOO [61I] are also given in the frames for comparison. 
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B. Direct detection 



111 this subsection we discuss the direct detection constraints on the effective scalar WIMP models described by 
Eqs. ((H) - (HI). Direct detection experiments measure the recoil energy of the atomic nuclei when the WIMPs elastically 
scatter off them. Because of the complex structure of nucleons we have to study the WIMP-nucleon interactions with 
another set of effective interactions in terms of the nucleon fields, by substituting the quark field operators in Eqs. ^ 
- ^ with the nucleon field operators TV (TV = p, n) and associating the coupling constants of these two sets of effective 
interactions by the form factors of nucleons. 

The velocity of the WIMP near the Earth is considered to be of the same order as the orbital velocity of the Sun, 
V ~ 0.001c, and we can safely calculate the scattering cross sections in the low velocity limit. In this limit, only 
the time component of (j>H^(f> and the spacial components of Nj^^j^N survive, thus the VA interaction leads to an 
elastically scattering cross section with velocity suppression. On the other hand, Nj^N vanishes in this limit so that 
the SP interaction is not constrained by direct detection either. Therefore there is no measurable spin-dependent 
(SD) WIMP-nucleon scattering. This situation can be also understood by realizing that scalar WIMPs have no spin 
structure. The scattering cross sections induced by the remaining S and V interactions are both spin-independent 

(SI). ^ ^ .... . . . . 

The calculation of scattering cross sections is much similar to our previous calculation for fermionic DM [50|, so 
we will just outline the result. To compare with the result of CDMS [s^ and XENON [6l|, the WIMP-nucleon cross 
sections are given as follows: 



2 



Scalar int. : .s^a^ = ^^^m'^^^,^)^ (^) ' (16) 

Vector mt. : ctv.^at = , , , 17 

where Fjv are the induced coupling constants of effective WIMP-nucleon interactions. Fn are related to the couplings 
to quarks Fq by form factors. For scalar interaction, 

Fs^N^ E Ps^X—+ E ^s./g" — . (18) 
^ nia ^ rria 

q—u,d,s ^ q—c.b,t 

where the nucleon form factors are fP = 0.020 ± 0.004, = 0.026 ± 0.005, /f = 0.118 ± 0.062, = 0.014 ± 0.003, 



= 0.036 ± 0.008, = 0.118 ± 0.062 [7^,1731, and fq ^ ^{l - fu - fd - /f ) for heavy quarks. For vector 



interaction, 

Pv,p = 2Fv,u + F\iM : Pv,n = Fv,u + 2i^v,d 7 (19) 

which reflect the valence quark numbers in the nucleons. For the case of real scalar WIMPs, the vector interaction 
vanishes, and the scattering cross section of scalar interaction should include an additional factor 4 comparing to the 
case of complex scalar WIMPs. 

In Fig. [31 we show the predicted elastic scattering cross sections between complex scalar WIMP and nucleo n by 
different kinds of effective interactions. Then given the experimental bounds set by CDMS II [s^ and XENON 100 [6l| . 
we can derive the upper bounds of coupling constants for different WIMP masses. The SP and VA interactions are 
not constrained by direct detection experiments, as stated above. Additionally, due to a difference of factor in 
Eqs. (|16|) and (|17|) . the predicted scattering cross section curves slope downward for scalar interactions, while they are 
nearly horizontal for vector interactions. This difference is a reflection of the different dimensions of the two effective 
operators. 



C. Indirect detection 



Compared with the direct search method for the WIMP DM at underground laboratories on Earth, an indirect 
detection method concerning astrophysical effects of DM is relatively independent and complementary. This method is 
used to look for the DM annihilation or decay products which include gamma rays, neutrinos, positrons, antiprotons 
and so on. These particles can be detected by satellite, ground and underground cosmic ray experiments. When 
propagating in the Galaxy, the charged products are deflected by Galactic magnetic fields and interact with interstellar 
medium. Therefore we have to consider the propagation process of the charged particles in order to compare predictions 
with observations. 
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FIG. 4. The 3a upper bounds on the coupHng constants Ff from the PAMELA p/p ratio [63 in each complex scalar WIMP 
effective model. In the left (right) column, results are given for the case of universal couplings {Ff cx m/). 



The propagation equation of cosmic rays in the Galaxy can be written down as [74 



d_ 

dp 



pip- 



|(V.V„ 



- (20) 

Tf Tr 



where "0 = Tp(r,p,t) is the number density of cosmic-ray particles per unit momentum interval, Q(r,p) is the source 
term, D^^ is the spatial diffusion coefficient, Vc is the convection velocity, Dpp is the diffusion coefficient in momentum 
space describing the reacceleration process, p = dp/dt is the momentum loss rate, Tf and Tr are the time scales for 
fragmentation and the radioactive decay, respectively. Following our previous work [501 consider the antiproton- 
to-proton flux ratio p/p measured by the satellite-borne experiment PAMELA [g^ to constrain the effective models. 
To proceed we solve the propagation equation 
to the PAMELA p/p data. 



20)) by the numerical code GALPROP [7^ and compare the solution 
The source term of p induced by the annihilation of complex scalar WIMPs is given by 



Qann (r, E) 



dNf 
dEe 



u, d, s, c, b, t, 



(21) 



where (crannw)tot total annihilation rate of WIMPs to quarks, Bq is the branching ratio of each qq channel, 

pir) is the mass density distribution of the Galactic DM halo, [dNp/ dEp)q is the number per unit energy interval of 
the antiprotons produced by a WIMP pair-annihilation in the qq channel. In particular, {dNp/dEp)q is calculated by 
the Monte Carlo code PYTHIA [t^. For real scalar WIMPs, the factor ((Tannw)tot in Eq. ^ is modified to 

(canni')tot /'^-^^4> siucc the two annihilating WIMP are indistinguishable. 

The NFW profile [t^ is taken to describe p{r) with the characteristic density ps = 0.334 GeV cm^ and the scale 
radius = 20 kpc, which leads to a local DM density /5(8.33 kpc) = 0.4 GeV cm^"^. The DM particles in the 
Galactic halo obey the Maxwell-Boltzmann velocity distribution f{vo) = {M^/2TTkBT)'^^^ exp{—M^VQ/2kBT). Their 

jal. In the calculation 



velocity dispersion vq 



y^3T/M^ is taken to be the canonical value 270 km s 
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of the p/p spectrum with GALPROP, the Galaxy propagation model with diffusion and convection is adopted and 
the half-height of the Galaxy propagation halo is set to be 4 kpc. The PAMELA p/p spectrum is best fitted by 
astrophysical background with a solar modulation potential of $ = 335 MV, which gives the minimal . Under 
such potential if the contribution of DM is added in, the resulted spectrum deviates from the background and gives 
a larger . Thus for the fixed in both the cases of universal couplings and Ff (x uif for each type of effective 
interactions, we can derive the 3cr upper bounds on ((Jannw)tot' which eventually set the bounds on effective coupling 
constants Fj, as shown in Fig. HI 

In Fig. m the upper bounds set for S and SP interactions and those for V and VA interactions are nearly identical 
respectively, with slight differences at ~ nit in the case oi Ff (x uif. The reason is the same as that explained in 
Subsection III Al It is also worth noting that the upper bound curves for the operators with different dimensions have 
different trends of slope, due to the additional factor Af^ in the annihilation cross sections of dimcnsion-6 operators 
(V and VA). For the case of Ff ex ti/, downward bends at A'/^ ^ mt in the curves appear again, because of the 
— > ti threshold effect. 
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FIG. 5. Combined constraints on coupling constants Fj of complex scalar WIMPs with scalar (S) interaction from relic density, 
direct detection experiments of CDMS II and XENON 100, PAMELA p/p ratio, and validity of effective theory. The yellow 
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FIG. 6. Combined constraints on coupling constants Ff of complex scalar WIMPs with vector (V) interaction from relic density, 
direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. 
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FIG. 7. Combined constraints on coupling constants Ff of complex scalar WIMPs with scalar-pseudoscalar (SP) interaction 
from relic density, PAMELA p/p ratio, and validity of effective theory. 
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FIG. 8. Combined constraints on coupling constants Ff of complex scalar WIMPs with vector-axialvector (VA) interaction 
from relic density, PAMELA p/p ratio, and validity of effective theory. 



D. Validity region of effective models and combined constraints 



In this subsection, we will discuss the validity region "where the method of effective theory can be used. Since our 
previous analysis is made in both the cases of universal couplings and Ff oc nif, we consider them case by case in the 
follo"wing. 

• In the case of universal couplings, for dimension-5 operators S and SP (dimension-6 operators V and VA), Ff 
can be written as Ff/\/2 = a/A {Ff/V2 = a/A^), where A is the cutoff energy scale and a is the coupling 
of the fundamental theory beyond A, which may be of order 1. The transfer momentum of the annihilation 
process cfxj)^ — > // must be well below the cutoff, i.e., 2Mcj, <C A, so that the effective theory can be used. On 
the other hand, a weakly coupled UV completion of the effective theory usually requires a < 4tt such that the 
perturbative method is valid [43,[73|- From the above 3 relations, we obtain 



for S and SP interactions, 
for V and VA interactions. 



(22) 
(23) 
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In the case of F/ cx nif, for dimension-5 operators S and SP (diniension-6 operators V and VA), Ff can be 
written as = amf/K^). Likewise, we still have 2M^ ^ A and a < 4tt. From the 

above 3 relations, we obtain 

Ff \f2-K 

— <C for S and SP interactions, (24) 

F( TT 

^ — -, for V and VA interactions. (25) 



TABLE I. A summary for complex scalar WIMPs with various effective interactions. The excluded regions of given by 
direct and indirect experiments are indicated. 





Universal coupling 




Interaction 


Direct detection 


PAMELA vIp 


S 


Excluded ~ 10 GeV - above 1 TeV Excluded ~ 10 - 73 GeV 


V 


Excluded M4, ~ 10 GeV - above 1 TeV 


Not sensitive 


SP 


Not sensitive 


Excluded ~ 10 - 73 GeV 


VA 


Not sensitive 


Not sensitive 




Ff oc mj 




Interaction 


Direct detection 


PAMELA p/p 


S 


Excluded ~ 10 - 168 GeV 


Excluded ~ 10 - 73 GeV 


V 


Excluded ~ 10 - 100 GeV 


Not sensitive 


SP 


Not sensitive 


Excluded ~ 10 - 73 GeV 


VA 


Not sensitive 


Not sensitive 



Considering the validity conditions p2p - ([25]) altogether with the other phenomcnological constraints, we derive 
the combined constraints on the effective models. In Figs. [5] and [5] we show the combined constraints on coupling 
constants Ff of complex scalar WIMPs with S and V interactions from relic density, direct detection experiments 
of CDMS II and XENONIOO, PAMELA p/p ratio and validity of effective theory, while in Figs. [7] and H we show 
the combined constraints on Ff with SP and VA interactions from relic density, PAMELA p/p ratio and validity of 
effective theory. The invalid parameter space of effective field theory is denoted by the yellow region in each figure. 
The constraints are much similar between S and SP interactions, and so are those between V and VA interactions, 
except that the S and V interactions are constrained by direct detection experiments, while the SP and VA interactions 
are not. 

There are some interesting features in Figs. [5] - [H In the case of universal coupling for S and V interactions, the 
WIMP with < 1 TeV arc strongly excluded by direct detection experiments, while in the case Ff oc nif, the 
direct detection exclude a much smaller range of M^. The indirect detection constraints also rule out the WIMP with 
< 70 GeV for S and SP interactions, while the constraints becomes much weaker for V and VA interactions. The 
reason is that for V and VA interactions the annihilation rate pH)) and ([T^ are of order 0{T/M^), which lead to looser 
bounds by indirect detection because of the much lower WIMP temperature today comparing with that of the WIMP 
freeze-out epoch. Additionally, in the case of oc m/ the effective interactions are valid only up to ~ TeV. A 
larger WIMP mass in these cases may implies violation to our previous assumptions, such as the existence of a light 
intermediate state, or the effect of resonances or coannihilations. 

If the effective coupling is much weaker than that derived from the observed relic density, the thermal DM production 
may overdose the Universe. Therefore the DM models in which the direct or indirect detection constraints are stronger 
than that of relic density should be excluded, or else some exotic entropy generation processes should occur after DM 
froze out. As a summary of the study on the scalar WIMP, in Tabic HI we indicate the excluded regions of given 
by direct and indirect experiments for complex scalar WIMPs with various effective interactions. 
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III. VECTOR DARK MATTER 

111 this section, we will discuss the case that DM consists of complex vector WIMPs (X^ and its antiparticle X'^*). 
Notable differences between this case and the case of real vector WIMPs will be remarked below. Following the 
procedure in which we discuss scalar WIMPs in Sectionllll we construct effective interaction operators oi X'^*-X'^-f-f 
type between the WIMPs and the SM fermions, which are limited only by the requirements of Hermiticity, Lorentz 
invariance and CPT iiivariance. 

The assumptions listed in Section |ll] are also adopted. Thus the effective interaction term of Lagrangian between 
two complex vector WIMPs {X^ and X^*) and two SM fermions (/ and /) is given by only one of the following 
expressions: 

Scalar int. (S) : = E ^K^'^ff, (26) 

Vector int. (V) : = ^ ^{X*Mxnh^f, (27) 

Tensor int. (T) : = ^i{X;X, - X,:X^)/a^V, (28) 

Scalar-Pseudoscalar int. (SP) : £sp = > ^ ^^^X*X^'fij5f, (29) 

Vector- Axialvector int. (VA) : £vA = > ' ^^^^{X;iB^X'')fY'j5f, (30) 

v2 





/ 


Ct = 






f 


^SP 








-CvA 


= E 






^V = 






f 




= E 











Alternative Vector int. (V) : = } ] ^^£'"''"'iX*KXp)f"f^f, (31) 

v2 



Alternative Vector- Axialvector int. (VA) : = ^ ^^^^e^'^P" {X*K Xp)f-f^-f5f, (32) 

v2 

Alternative Tensor int. (T) : ^f = Yl ^X^X^ ~ X;Xp)fapaf, (33) 

where the sum of / is over all the SM fermions, and the effective coupling constants Kf are real numbers. For 
the dimension-5 operators (S, SP, T and T), the coupling constants Kf have mass dimension of —1, while for the 

dimension-6 operators (V, VA, V and VA), Kf have mass dimension of —2. Operators with more derivatives on 
the fields are not considered, since such operators arc more suppressed at low energy scale due to their higher mass 
dimensions. Note that for the case of real vector WIMPs, the V, T, VA and T interactions vanish. In the following, 
we will discuss the phenomenological constraints on Kf given by the recent results of the DM relic density, direct and 
indirect detection experiments. 

A. Vector WIMP annihilation and relic density 

In order to calculate relic density and DM contribution in cosmic ray, the annihilation cross sections of WIMP- 
antiWIMP to SM fermion-antifermion pairs are needed. For each effective model listed above, the result is given 

by 



1 ^f^^sjV /5-4m2 (s-4m2)(s2-4A4s-f 12A4) 



f \ V ^ / V ^^"^ 

/Kv fV I (s + 2m'^f){s^ ~ AMls + 12Mi) 



4327rA/4 ^ V V2 



1 sr^fKTfV s-imj 
'^^-^ - 216;^ ^ j '^V^^4Mf 



+ {8m} + 4:M^)s - l&m)Ml - 8A4 - 160 ^ ^ 



(36) 
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FIG. 9. The predicted coxipling constants Kf as functions of the vector WIMP mass Mx, fixed by the observed relic density, 
noMh^ = 0.1109 ±0.0056 0]. Frames in the left (right) column are results given for the case of universal couplings [Kf cx m/). 



CSP, 



CVA, 



2887rM 



V2 



c/1 



s - 4M| 



1 



4327rA/_y 



VA,/ 

^/2 



CfJis-4m}){s^4Ml) 



1 



V, ann IQgj^M 



VA.ann IOSttM 



1 



{s - UlDis + 2mj) 



s - 4(Mi + mj) + 28 



rnjMl 



V2 



I s — 4mj 
s - 4Ml 



+ 4(M|- - mj)s + 32mjMl - SM^ + 128 



(37) 
(38) 

(39) 

(40) 

(41) 



where Mx is the WIMP mass. For the case of real vector WIMPs, each non-vanishing annihilation cross section 
should include an additional factor of 4 because the WIMP and its antipartncr arc identical. After taking the thermal 
average, we obtain the annihilation rates, i.e., 



^">^T^E 



12.^V^/2 j '^V Mi 



1 - 



(-V.a„n.)c.^E(^ 



Ml) 


•5(- 


mi 
h 7 — — 
Ml 









(42) 
(43) 
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FIG. 10. The predicted thermal rehc density (dashed hnes) of complex vector WIMPs with S, V, T, T, V and VA interactions, 
respectively. In the upper two frames, results are shown for the case of Kf oc m/ with S and V interactions. In the remaining 
frames, results are given for the case of universal couplings. The narrow horizontal solid band shows the range of the observed 
DM relic density, ^dm/i^ = 0.1109 ± 0.0056 
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Note that in the case of V interaction, {(Jannv) vanishes in the zeroth and the first orders of T/Mx, thus we retain 
the term of order T'^/M^. 

The procedure to obtain the rehc density of vector WIMPs is similar to that in the case of scalar WIMPs. A 
notable difference is that the degree of freedom for a vector WIMP is g = 3. Besides, since (aann'^) of V interaction 
(|47| is of order T'^/M'^, Eqs. (|14p and (|15p need to be extended. If CTann'^ in the lab frame is expanded to be 
(TannW — ^0 + ^1^^ + o.2v'^ + 0{v^), wc have ((Tannw) ~ oq + 6aix~^ + (60a2 — 27ai)x~^ + 0{x~^) with X = Mx/T. 
Then the factor (oq + SaixJ^) in the denominator of the right-hand side of Eq. should be modified to be 

[ao + SaixJ^ + (20a2 — 9ai)xj'^], and the factor (oq + 6aixJ^) in the square bracket of the right-hand side of Eq. ([15]) 

should be modified to be [ao -I- 6aixJ^ + (60a2 — 27ai)xJ^]. For the case of real vector WIMPs, the relic density just 
one half of that for the case of complex vector WIMPs given the same annihilation cross section. 

In Fig. [SI the predicted coupling constants Kf fixed by the WMAP observed relic density are shown as functions 
of Mx- The upper two frames are plotted for dimension-5 operators (S, SP, T and T), while the remaining two are 
plotted for dimension-6 operators (V, VA, V and VA). Two cases of coupling constants are considered again. In the 
left (right) column of the frames, results are given for universal couplings (Kf cx rrif). There are some interesting 
features in Fig. [HI 

• The curves for the case of Kf oc mj have sharp downward bends at Mx ^ rrit- The reason is the same as that 
in the case of scalar WIMPs. Even for the case of universal couplings, the curve of VA interaction also exhibits 
the same bending behavior because of <(t,tx v) cx rni. 

^ \ VA, ann / J 

• There are two pairs of nearly identical curves, S ~ SP and V ~ VA. This is because of the similarities of the 
annihilation rates corresponding to these two pairs of interactions, respectively. The leading terms of these two 
pairs of the annihilation rates tend to be the same when Mx is much larger than that of the final state fermions. 

• The curves of V and VA interactions are obviously lower than the other two curves of V and VA interactions. 
That is due to the forms of the annihilation rates. The leading terms of the annihilation rates for V and VA 
interactions are of order T/Mx, while that for V interaction is of order T'^/M^. The leading term of the 
annihilation rate for VA interaction is proportional to mjT/Mx, which is suppressed by the SM fermion masses 

771/ that are smaller than Mx almost in the whole region. On the other hand, the curves of T and T interactions 
are lower than those of S and SP interactions due to the larger zeroth order terms of the annihilation rates for 
T and T interactions. 

In Fig. [TUl we show the predicted thermal relic density JIdm^^ of complex vector WIMPs for S, V, T, T, V and VA 
interactions. The curves for SP and VA interactions are very similar to those for S and V interactions respectively, as 
discussed above. In the upper two frames, results are shown for the case of Kf oc m/ with S and V interactions. The 
remaining frames give the results for the case of universal couplings. Similar to the case of scalar WIMPs, the curves 
for dinicnsion-5 operators (S, SP, T and T) arc nearly horizontal in almost the whole mass range, while those for 
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dimension-6 operators slope downward (except for those of VA interaction, due to (cy^ ann^) ^ "^Z^" Additionally, 
if resonances, coannihilations, or annihilations to final states other than fermion-antifermion pairs are significant, the 
actual curves in Figs. and [TUl will be significantly lower than they are shown there. 



B. Direct detection and indirect detection 



In this subsection we discuss the direct and indirect detection constraints on the effective vector WIMP models 
described by Eqs. (^5]) ~ (|33p . In the low velocity limit, only the scattering cross sections of S, V, T and VA interactions 
survive, and others are suppressed by the low WIMP velocity. The S and V interactions between WIMPs and nuclei 
are spin- independent (SI), while the T and VA interactions are spin-dependent (SD). In order to be compared with 
the results of CDMS [H^fi^l and XENON [61|[62|, the WIMP-nuclcon cross sections arc given as follows: 



Scalar int. : era xn ~ ^ ( — ^ ) , (50) 



Vector mt.: ay.xjv = — =- ) , (51) 

TT [Mx + mN) \ V ^ 



lensor mt. : ctt xn = ^ — 1=- , (52 

TT {Mx + mN f \ V2 J 

VAint.: a^^,,= V ^ (53) 

^^^^^ n{Mx+mNf \ V2 J ' 

where {N = p, n) are the induced coupling constants of effective WIMP-nucleon interactions, related to the 
couplings of WIMP-quark Kq by form factors. For S and V interactions, the relations between Kn and Kq are the 
same as those between Fn and Fq (Eqs. (fT8|) and ([T9| ). respectively. For VA and T interactions, 

i^VA.Ar= ^VA.gAf, (54) 

q—u.d,s 

Kt,n= ^T,,Af, (55) 

q—u.d,s 

with the form factors 0.842 ± 0.012, A^ = -0.427 ± 0.013, Af = -0.085 ± 0.018 [zl, A^J = A^, A^ A^, 

A" = Af , which reflect the contributions of quark components to the nucleon spin. Note again for real vector WIMPs, 

the V and T interactions vanish, and the scattering cross sections of S and VA interactions should include an additional 
factor 4 comparing to the case of complex vector WIMPs. 

To compare with the searching bounds of direct detection experiments, we plot the predicted cfxn vs. Mx curves 
with fixed couplings. In Fig. [11] the spin-dependent WIMP-proton cross sections (dashed lines) for S and V interactions 
are shown. In the left (right) column, results are given for universal couplings {Kf oc to/). The upper bounds set by 
CDMS II [iil and XENONIOO ^^aie also given in the frames. In Fig. [HJ the spin-dependent (SD) WIMP-neutron 
cross sections (dashed lines) for VA and T interactions are shown. In the left (right) column, results are given for 
universal couplings {Kf oc to/). The upper bounds set by CDMS [6^ and XENONIO [i^ are also given in the frames. 

From Figs. [TT] and I12[ we find that the experimental constraints are much stronger for SI interactions than those 
for SD interactions, since the scattering amplitudes between WIMPs and nucleons in the nuclei add coherently in the 
case of SI interactions. Additionally, the curves for dimension-5 operators (S and T) slope downward, while those for 
dimension-6 operators (V and VA) remain horizontal in most range of mass. This is due to an additional factor of 
M^ in Eqs. dH]) and comparing with Eqs. and (|5^ . 

The calculation of indirect detection constraints of vector WIMPs is totally similar to that of scalar WIMPs in 
Subsection III CI In Fig. [T31 we show the 3tT upper bounds on the coupling constants Ff from the PAMELA p/p 
ratio j67| in each effective complex vector WIMP model for the cases of universal couplings and Kf oc m/. 

Duo to the similar reasons, the bounds in Fig. 1131 have similar properties as the relic density bounds in Subsection 
nil Al (a) There are sharp downward bends in the curves at Mx ^ mt when (cTannw) oc or oc m^. (b) There are 
two pairs of interactions that have almost the identical bounds, S ~ SP and V ~ VA. (c) The curves of S and SP 
interactions lie well above those of T and T interactions, while the curves of V and VA interactions lie well below those 
of V and VA interactions. Again, we can find the explanations to these features in the annihilation rates Eqs. (|42p - 
Eq. (dHD. 
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FIG. 11. The spin-independent (SI) WIMP-proton cross sections (dashed hnes) for complex vector WIMPs with S and V 
interactions. In the left (rig ht) column, results are given for universal couplings (A'/ oc m/). The upper bounds set by CDMS 
II [S^ and XENONIOO [6l| are also given in the frames for comparison. 
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FIG. 12. The spin-dependent (SD) WIMP-neutron cross sections (dashed lines) for complex vector WIMPs with VA and 
T interactions. In the left (rig ht) column, results are given for universal couplings {Kf oc m/). The upper bounds set by 
CDMS [60| and XENONIO [62| are also given in the frames for comparison. 
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C. Validity region of effective models and combined constraints 

Since we arc using effective theory, we can only carry out our analysis below an energy cutoff scale A, which leads to 
validity bounds on the effective couplings Kf if the perturbative theory is assumed for UV completion, as we discussed 
in Subsection III Dl Then for the case of universal couplings, we have 

for S, SP, T and T interactions, (56) 

for V, VA, V and VA interactions. (57) 

On the other hand, for the case of Kf (x nif we obtain 

for S, SP, T and T interactions, (58) 

for V, VA, V and VA interactions. (59) 

Considering the validity conditions ([55)) - ([5^ altogether with other phenomenological constraints, we obtain the 
combined constraints on the effective vector WIMP models. In Figs. [U]- [121 the combined constraints of the effective 
models of S, V, T, V, VA and T interactions are shown, with the yellow regions denoting the invalid parameter spaces 
of effective field theory. The constraints on SP and VA interactions are very similar to those on S and V interactions 
respectively, except that the SP and VA interactions are not constrained by direct detection experiments. S and V 
interactions with universal couplings are stringently constrained by direct detection experiments, while the constraints 
for the case of Kf oc nif are much looser. The direct detection bounds for SD interactions (T and VA) is so much looser 
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than those for SI interactions (S and V) that they only exclude a region of small mass in the case of VA interaction 
with universal couplings. The constraints by PAMELA p/p ratio exclude the mass regions of Mx ^ 70 GeV for 
dimension-5 operators (S, SP, T and T), while those for dimension-6 operators (V, VA, V and VA) are much weaker. 
The reason is that for dimension-6 operators, the annihilation rates Eqs. (|33|), (|46p . (|Tf)) and (|3S|) are of order T/Mx 
or T^/M|-, which lead to much looser bounds by indirect detection because the WIMP temperature today is much 

lower than that of the WIMP freeze-out epoch. In the cases of V and VA interactions with Kf oc m/, large parts of 
the curves set by relic density lie in the invalid regions of effective theory, which may indicate effects of resonances, 
coannihilations, light intermediate states, or final states other than SM fermion pairs. 
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FIG. 14. Combined constraints on coupling constants Kf of complex vector WIMPs with scalar (S) interaction from relic 
density, direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. The 
yellow region denotes the invalid parameter space of effective field theory. The left (right) frame is shown for the case of universal 
couplings {Kf ocmf). The constraints from relic density, p/p ratio and validity of effective theory for scalar-pseudoscalar (SP) 
interaction are very similar to those for scalar interaction, but direct detection experiments are not sensitive to SP interaction. 




FIG. 15. Combined constraints on coupling constants Kf of complex vector WIMPs with vector (V) interaction from relic 
density, direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. The 
constraints from relic density, p/p ratio and validity of effective theory for vector-axialvector (VA) interaction are very similar 
to those for vector interaction, but direct detection experiments are not sensitive to VA interaction. 

As a summary of the study on the vector WIMP, in Table [TTl the excluded regions of Mx given by direct and 
indirect experiments for complex vector WIMPs with various effective interactions are shown. 



IV. SPIN-3/2 DARK MATTER 

In this section, we will extend our discuss to the case that DM consists of spin-3/2 fcrmionic WIMPs. The spin-3/2 
WIVIP can be described by a vcctor-spinor Xa which carries both a vector index fi and a spinor index a. Then a free 
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FIG. 16. Combined constraints on coupling constants Kj of complex vector WIMPs with tensor (T) interaction from relic 
density, direct detection experiments of CDMS and XENONIO, PAMELA p/p ratio, and validity of effective theory. 
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FIG. 17. Combined constraints on coupling constants of complex vector WIMPs with alternative vector (V) interaction 
from rehc density, PAMELA p/p ratio, and validity of effective theory. 
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FIG. 18. Combined constraints on coupling constants Kj of complex vector WIMPs with alternative vector-axialvector (VA) 
interaction from relic density, direct detection experiments of CDMS and XENONIO, PAMELA p/p ratio, and validity of 
effective theory. 



WIMP satisfies the rank-1 Rarita-Schwinger equations [7^| 

(i^-Mjx^ = 0, 7mX^ = 0, (60) 
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FIG. 19. Combined constraints on coupling constants Kf of complex vector WIMPs with alternative tensor (T) interaction 
from relic density, PAMELA p/p ratio, and validity of effective theory. 



TABLE II. A summary for complex vector WIMPs with various effective interactions. The excluded regions of Mx given by 
direct and indirect experiments are indicated. 





Universal coupling 




Interaction 


Direct detection 


PAMELA p/p 


S 


Excluded Mx ^ 10 GeV - above 1 TeV Excluded Mx ~ 10 - 74 GeV 


V 


Excluded Mx ^ 10 GeV - above 1 TeV 


Not sensitive 


T 


Not sensitive 


Excluded Mx ~ 10 - 71 GeV 


SP 


Not sensitive 


Excluded Mx ~ 10 - 73 GeV 


VA 


Not sensitive 


Not sensitive 


V 


Not sensitive 


Not sensitive 


VA 


Excluded Mx ~ 10 - 170 GeV 


Not sensitive 


f 


Not sensitive 


Excluded Mx ~ 10 - 72 GeV 




Kf (X ruf 




Interaction 


Direct detection 


PAMELA p/p 


S 


Excluded Mx ~ 10 - 172 GeV 


Excluded Mx ~ 10 - 74 GeV 


V 


Excluded Mx ~ 10 - 157 GeV 


Not sensitive 


T 


Not sensitive 


Excluded Mx ~ 10 - 70 GeV 


SP 


Not sensitive 


Excluded Mx ~ 10 - 74 GeV 


VA 


Not sensitive 


Not sensitive 


V 


Not sensitive 


Not sensitive 


VA 


Not sensitive 


Not sensitive 


f 


Not sensitive 


Excluded Mx ~ 10 - 72 GeV 



where the spinor indices have been suppressed. From Eqs. (|60|) . it can be seen that each component of acts as 
a Dirac spinor, while all the four components together are constrained by the subsidiary condition 7^^%^ = 0. By 
introducing an orthogonal spin projection operator, the spin sum relations in momentum space can be derived [80ll8l|. 
i.e., 



J2 <ij>Ws{p) = it + AIx) [P"" Ip'^'P'^IpI.) , (61) 
J2 v^ipKiP) -kf- M^) (P"" - Ip^'P'^IpI.) , (62) 



where s is the spin index, and u^(p) (vsip)) the positive (negative) energy solutions of Eqs. (|60p . and P^'' 
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Although spin-3/2 field theories suffer from a common problem of non-renormalizability, it is safe to consider effective 
field theory below an appropriate energy cutoff. Following the assumptions in Sections and lllli we construct the 
effective interaction operators of yf-X^-f-f type between the WIMPs and the SM fermions as follows, 



Scalar int. (S) : = ^ ^^xjl, (63) 

Pseudoscalar int. (P) : Cp = J2 ^x''l5xjl5f, (64) 



V2 



Vector int. (V) : = ^ -^x'^XpHf, (65) 



Axialvector int. (A) : = E -^x'^l^XpHl^f, (66) 

Tensor int. 1 (Tl) : ^^i = ^ ^^X^^^^Xm^p-/, (67) 



Tensor int. 2 (T2) : ^^2 = ^|^z(x^x" - rx^)/>M-/, (68) 
Tensor int. 3 (T3) : -Cts = ^ ^|^£^'^'"^*(x;.X. - X.Xf.)f'^p.f, (69) 



Tensor int. 4 (T4) : Ct4^Y1 ^^^''"'"x^'^^'^Xr.f'Jpa.f, (70) 
Tensor int. 5 (T5) : Ct^^YI %^(x^75X'' " ri^X^f'^.^f. (71) 



/ 



V2 

Gt6. 



Tensor int. 6 (T6) : Cts ^WX. " X.l5X^)f^paf, (72) 



^/2 

Gsp,/_ 



Scalar-Pseudoscalar int. (SP) : £sp = X! -^^^X^Xm/Hs/. (73) 
Pseudoscalar-Scalar int. (PS) : Cps=Y ^y^X^'iTsX^//- (74) 



/ 



Vector-Axialvector int. (VA) : CvA = J2 ^^X^t'^Xm/t-^Ts/, (75) 

Axialvector- Vector int. (AV) : -^^X^'l'lsxjl.f, (76) 

where the sum of / is over all the SM fermions, and the effective coupling constants G/ are real numbers with mass 
dimension of —2. Since each component of x'' looks like a Dirac spinor, the S, P, V, A, Tl, T4, SP, PS, VA and AV 
interactions are analogous to the S, P, V, A, T, T, SP, PS, VA and AV interactions defined in our previous work on 
Dirac fermionic WIMPs [H^l, respectively. Even though the spin-3/2 fermionic WIMP has more degrees of freedom, 
it is not important in the calculation of phenomenological constraints. Thus we can expect the phenomenological 
constraints of spin-3/2 WIMPs are also much similar to those of Dirac fermionic WIMPs. For this reason, we just 
give a brief discussion on spin-3/2 WIMPs in the following. 

The annihilation cross sections of WIMP-antiWIMP to SM fermion-antifermion pairs are given by 
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Taking the thermal average, we obtain 
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When calculating the relic density of spin-3/2 WIMPs, it is worth noting that the degree of freedom for a spin-3/2 
fermionic WIMP is g = 4. 

In the low velocity limit, the scattering cross sections of P, T3, T4, T5, T6, SP, PS, VA and AV interactions between 
WIMPs and nuclei vanish. In the remaining interactions, the S and V interactions are spin-independent, while the A, 
Tl and T2 interactions are spin-dependent. The WIMP-nucleon cross sections are given by 



Scalar int. 
Vector int. 



tt{M^ + hinY V \/2 
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mlMl 



G 



V,N 
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(107) 
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5m2,M2 fGA,NV 



Axialvector int. : cta yw = „ , . , — ttt — ^ , (108) 



Tensor int. 1 : (Tti vat = , " ^ — ^ , (109) 

Tensor int. 2 : axz xiv = , ^ „ , (110) 

where Gat (A^ = p, n) are the induced coupHngs of effective WIMP-nucleon interactions, related to the couphngs of 
WIMP-quark Gq by form factors. For S, V, A and Tl (or T2) interactions, the relations between Gat and Gq are 
similar to Eqs. (HH), (HI]), ([Ml) and respectively. Since ah the spin-3/2 WIMP effective models 1^ - have 
effective couplings G/ with mass dimension of —2, their validity regions are similar to those indicated in Eqs. ([23]) 
and (123). 




FIG. 20. Combined constraints on coupling constants G/ of spin-3/2 WIMPs with scalar (S) interaction from relic density 
direct detection experiments of CDMS II [53| and XENONIOO ^ij], PAMELA p/p ratio [631, ^-^d validity of effective theory. 
The yellow region denotes the invalid parameter space of effective field theory. The left (right) frame is shown for the case 
of universal couplings (G/ oc m/). The constraints from relic density, p/p ratio and validity of effective theory for scalar- 
pseudoscalar (SP) interaction are very similar to those for scalar interaction, but direct detection experiments are not sensitive 
to SP interaction. 




FIG. 21. Combined constraints on coupling constants G/ of spin-3/2 WIMPs with pseudoscalar (P) interaction from relic 
density, PAMELA p/p ratio, and validity of effective theory. The constraints from relic density, p/p ratio and validity of 
effective theory for pseudoscalar-scalar (PS) interaction are very similar to those for pseudoscalar interaction. 
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FIG. 22. Combined constraints on coupling constants G/ of spin-3/2 WIMPs with vector (V) interaction from relic density, 
direct detection experiments of CDMS II and XENONIOO, PAMELA p/p ratio, and validity of effective theory. The constraints 
from relic density, p/p ratio and validity of effective theory for vector-axialvector (VA) interaction are very similar to those for 
vector interaction, but direct detection experiments are not sensitive to VA interaction. 
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FIG. 23. Combined constraints on coup ling constants Gf of spin-3/2 WIMPs with axialvector (A) interaction from relic density, 
direct detection experiments of CDMS [6y| and XENONIO [63, PAMELA p/p ratio, and validity of effective theory. 
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FIG. 24. Combined constraints on coupling constants Gj of spin-3/2 WIMPs with tensor interaction 1 (Tl) from relic density, 
direct detection experiments of CDMS and XENONIO, PAMELA p/p ratio, and validity of effective theory. The constraints 
from relic density, p/p ratio and validity of effective theory for tensor interaction 5 (T5) are very similar to those for Tl 
interaction, but direct detection experiments are not sensitive to T5 interaction. 
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FIG. 25. Combined constraints on coupling constants Gj of spin-3/2 WIMPs with tensor interaction 2 (T2) from relic density, 
direct detection experiments of CDMS and XENONIO, PAMELA p/p ratio, and validity of effective theory. 




FIG. 26. Combined constraints on coupling constants Gj of spin-3/2 WIMPs with tensor interaction 3 (T3) from relic density, 
PAMELA p/p ratio, and validity of effective theory. 




FIG. 27. Combined constraints on coupling constants G/ of spin-3/2 WIMPs with tensor interaction 4 (T4) from relic density, 
PAMELA p/p ratio, and validity of effective theory. The constraints from relic density, p/p ratio and validity of effective theory 
for tensor interaction 6 (T6) are very similar to those for T4 interaction. 



The combined constraints on coupling constants Gf of the spin-3/2 WIMP effective models of S, P, V, A, Tl, T2, 
T3, T4 and AV interactions from relic density d], direct detection experiments of CDMS and XENON [Ulel], 

and PAMELA p/p ratio l67[ are shown in Figs.HOl"!!!! with the yellow regions denoting the invalid parameter spaces 
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FIG. 28. Combined constraints on coupling constants G/ of spin-3/2 WIMPs with axialvector-vector (AV) interaction from 
relic density, PAMELA p/p ratio, and validity of effective theory. 



of effective field theory. The constraints on SP, PS, T5, T6 and VA interactions are very similar to those on S, P, 
Tl, T4 and V interactions respectively, except that the SP, T5 and VA interactions are not constrained by direct 
detection experiments. As stated above, the constraints on S, P, V, A, Tl, T4, SP, PS, VA and AV interactions are 
much similar to those on the analogous interactions of Dirac fermionic WIMPs in our previous work [soj , respectively. 
In this work, however, we adopt the latest constraint of XENONIOO SI direct detection [6l|, which is stronger than 
that adopted in our previous work [53|. Thus the constraints on S and V interactions of spin-3/2 WIMPs are more 
stringent than those for Dirac fermionic WIMPs we calculated before. On the other hand, the T2, T3, T5 and T6 
interactions have no analogue of Dirac fermionic WIMPs. From Figs. [25l and [26l we can see that the direct detection 
constraints on T2 and T3 interactions are so loose that only a little region of is excluded by SD direct detection 
in the case of T2 interaction with universal couplings. In addition, the constraints by PAMELA p/p ratio exclude the 
mass regions of < 70 GeV for T5 and T6 interactions, while it is not sensitive to T2 and T3 interactions because 
their annihilation rates ([QT]) and (|98)) are both of order T /M^. As a summary of the study on the spin-3/2 WIMP, 
in Table Uni the excluded regions of given by direct and indirect detection experiments for spin-3/2 fermionic 
WIMPs with various effective interactions are shown. 



V. CONCLUSIONS 



In this work, we give a general analysis of the 4-particle interaction between SM particles and DM which consists 
of scalar, vector or spin-3/2 WIMPs. The most general forms of the 4-particle operators up to dimension 6 have 
been considered. We find that for scalar, vector and spin-3/2 DM, the constraints from DM relic density, DM direct 
and indirect detection are complementary to each other. Thus the comparison among different kinds of experimental 
results gives us a complete picture about the current DM searches. In general, the constraints from SI direct detection 
are the most stringent, while those from SD direct detection are quite weak. On the other hand, for light DM (whose 
mass < 70 GeV) the cosmic-ray p/p data can be more sensitive than direct detection or relic density in some cases. 

Assuming one operator dominates the effective interaction between DM and SM fermions, we find that in some 
cases the constraints are so strong that the Universe will be overdosed by DM thermal production. If the standard 
cosmology is still retained, the DM models in such cases should be excluded, which are indicated in Tables HI HIl and Hill 
In the case of scalar DM, recent direct detection experiments exclude some regions for S and V interactions, while 
the PAMELA p/p ratio excludes some small regions (< 70 GeV) for S and SP interactions. In the case of vector 

DM, recent direct detection experiments exclude some Mx regions for S and V interactions, and for VA interaction 
only with universal couplings. The PAMELA p/p ratio, however, excludes some small Mx regions (< 70 GeV) for 
S, T, SP and T interactions, and most of these interactions cannot be excluded by direct detection. In the case of 
spin-3/2 DM, the constraints on most interactions (S, P, V, A, Tl, T4, SP, PS, VA and AV) are much similar to 
those on their analogous interactions of Dirac fermionic DM in our previous work [s^. Besides, the constraints on 
T2 and T3 interactions are so weak that only direct detection can exclude a little region of in the case of T2 
interaction with universal couplings. For T5 and T6 interactions, some small regions (< 70 GeV) are excluded 
by the PAMELA p/p ratio. Among the scalar, vector and spin-3/2 DM effective models, there are still some effective 
interactions to which the recent DM direct and indirect search experiments are not sensitive at all. 
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TABLE III. A summary for spin-3/2 WIMPs with various effective interactions. The excluded regions of M-^ given by direct 
and indirect experiments are indicated. 
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